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We consider an "impurity" with a spin degree of freedom coupled to a finite reservoir of non-interacting elec- 
trons, a system which may be realized by either a true impurity in a metallic nano-particle or a small quantum 
dot coupled to a large one. We show how the physics of such a spin impurity is revealed in the many-body 
spectrum of the entire finite-size system; in particular, the evolution of the spectrum with the strength of the 
impurity-reservoir coupling reflects the fundamental many-body correlations present. Explicit calculation in 
the strong and weak coupling limits shows that the spectrum and its evolution are sensitive to the nature of the 
impurity and the parity of electrons in the reservoir. The effect of the finite size spectrum on two experimental 
observables is considered. First, we propose an experimental setup in which the spectrum may be conveniently 
measured using tunneling spectroscopy. A rate equation calculation of the differential conductance suggests 
how the many-body spectral features may be observed. Second, the finite-temperature magnetic susceptibil- 
ity is presented, both the impurity susceptibility and the local susceptibility. Extensive quantum Monte-Carlo 
calculations show that the local susceptibility deviates from its bulk scaling form. Nevertheless, for special 
assumptions about the reservoir - the "clean Kondo box" model - we demonstrate that finite-size scaling is 
recovered. Explicit numerical evaluations of these scaling functions are given, both for even and odd parity and 
for the canonical and grand-canonical ensembles. 

PACS numbers: 73.23.Hk, 73.21.La, 72.10.Fk 



I. INTRODUCTION 

The Kondo problem describes a single magnetic impurity 
interacting with a sea of electrons. 1 At temperatures T of 
the order of or less than a characteristic scale, Tk, the dy- 
namics of the impurity and the sea of electrons become in- 
extricably entangled, thus making Kondo physics one of the 
simplest realizations of a strongly correlated quantum sys- 
tem. In its original context, the impurity was typically an 
element of the 3d or 4/ series of the periodic table, embed- 
ded in the bulk of a metal such as Cu with s conduction elec- 
trons. With the subsequent development of fabrication and 
control of micro- and nano-structures, it was pointed out 2,3 
that a small quantum dot with an odd number of electrons - 
small enough that its mean level spacing As is much larger 
than the temperature - could be placed in a regime such that 
it behaves as a magnetic impurity. 4 ' 5 '* The first experimen- 
tal implementations of this idea were naturally made by con- 
necting the "magnetic impurity" formed in this way to macro- 
scopic leads i 6 ' 7 ' 8 - 9 The flexibility provided by the patterning 
of two dimensional electron gas makes it possible, however, 
to design more exotic systems, by connecting the small mag- 
netic impurity dot to larger dots playing the role of the elec- 
tron reservoirs. Schemes to observe, for instance, two-channel 
SU(2)i2>iI or SU(4) 12 i 13 i 14 ' 15 i 16 i 17 i 18 i 19 Kondo have been imple- 
mented. 

When the bulk electron reservoir of the original Kondo 
problem is replaced by a finite reservoir, two energy scales 
are introduced: the Thouless energy associated with the 
inverse of the time of flight across the structure, and the mean 
level spacing Ar i 20 - 21 - 22 A natural question which arises is 
therefore how these two new scales affect the Kondo physics 



under investigation. 

Because a quantum impurity problem has point-like inter- 
actions, the local density of states Pioc(e) = X)a |0a (0) | 2 <5(e — 
e Q ) completely characterizes the non-interacting sea of elec- 
trons (e a and (f> a are the one-body eigenvalues and eigenfunc- 
tions of the reservoir). For T, Xk 3> -&rh, Ar, thermal smear- 
ing washes out the effects of both mesoscopic fluctuations 
and the discreteness of the reservoir spectrum. Indeed, in this 
regime, one may safely approximate p\ oc by a constant po\ the 
impurity behaves in much the same way as if it were in an infi- 
nite reservoir. In contrast, when X, Xk < -Erh, Ar, the impu- 
rity senses the finiteness of the reservoir through the structure 
of pioc(e). The presence of these new energy scales (which 
are ubiquitous 2021 in reservoirs made from quantum dots) is 
hence an essential and interesting part of Kondo physics in 
nano-systems and deserves to be understood thoroughly. 

The implications of a finite Thouless energy, and of the as- 
sociated mesoscopic fluctuations taking place in the energy 
range [A,£/rh], nave t> een investigated mainly in the high 
temperature range X;^ Xk, where a perturbative renormaliza- 
tion group approach is applicable 23 ,24 ' 25 ' 26 - 27 (see also related 
wor k 28 ' 29 ' 30 in the context of weakly disordered system). Less 
is known about the implications of mesoscopic fluctuations in 
the temperature range X < Tk. 

There is on the other hand already a much larger body of 
work concerning the "clean Kondo box" problem, 31 namely 
the situation where mesoscopic fluctuations are ignored (or 
absent as may be the case in some one dimensional mod- 
els), and only the existence of a finite mean level spacing is 
taken into account. Simon and Affleck 32 and Cornaglia and 
Balseirc 33 have, for instance, considered how transport prop- 
erties are modified if one dimensional wires of finite length 
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FIG. 1: (Color online) A double dot system, coupled very weakly 
to leads (LI, L2). In the Coulomb blockade regime, the small dot 
S behaves like a magnetic impurity that is coupled to a finite reser- 
voir R provided by the large dot. The leads are used to measure the 
excitation spectrum of the system. 

are inserted between the macroscopic leads and the quantum 
impurity. Ring geometries ^ 34 ' 35,36 ' 37 including the configura- 
tion corresponding to a two channel Kondo effect, 37 have also 
been investigated. 

The basic Kondo box configuration, namely a quantum im- 
purity connected to an electron reservoir with a finite mean 
level spacing, turns out to be already a non-trivial problem 
and so has been investigated by various numerically intensive 
techniques such as the non-crossing approximation 31 or the 
numerical renormalization group. 38 In Refs.HH and|37| it was 
pointed out, however, that as only the regime T <C is af- 
fected by the finiteness of A, a lot of physical insight could be 
obtained by the analysis of the low energy many-body spec- 
trum of the Kondo box system (i.e. the ground state and first 
few excited states). An analysis of this low energy many- 
body spectra and of an experimental setup in which it could 
be probed was given in Ref.[39l 

In this article, we would like on the one hand to provide 
a more detailed account of some of the analysis sketched in 
Ref.HH and furthermore to present an additional physical ap- 
plication, namely the low temperature magnetic response of 
the Kondo box system. See also Ref. 40 for an analysis of the 
addition energy of a Kondo box. 

Since our focus is the consequences of a finite Ar, we con- 
sider the simplest possible configuration: a double dot system 
with a small dot acting as the magnetic impurity and a larger 
one playing the role of the electron reservoir, as illustrated in 
Fig.Q] The Hamiltonian describing this double-dot system is 

H R . S = ^cLc aa + E C (N R - nf) 2 + H*i A . (1) 

Ota 

Here cj^ creates a state </> Q (r) with spin a and energy e Q 
which is an exact one-body energy level in the bigger quantum 
dot R. These states include all the effects of static disorder and 
boundary scattering. N R is the number operator for electrons 
in the reservoir, the dimensionless gate voltage applied to 
the large dot R, and Ec its charging energy. As the charg- 
ing energy is the leading interaction for electrons in a finite 
system, we shall neglect all other interactions among the elec- 



trons on R. (See, e.g., Ref. 40 4l1l42l for work that includes 
interactions among electrons in R.) The last term in Eq. (Q} 
contains the description of the small dot and the interaction 
between the dots. 

We consider two models for the magnetic impurity quan- 
tum dot and its interaction with the reservoir R. For most of 
this paper, we use a "Kondo"-like model, which therefore in- 
cludes charge fluctuations only implicitly. In this case, the 
smaller quantum dot is represented by a spin operator S. The 
interaction with the screening reservoir R is given by the usual 
Kondo interaction, 

if£ t — J S ■ s(0) , (2) 

describing the anti-ferromagnetic exchange interaction be- 
tween the two dots, with s(0) — \.fl a ~cr aa' foa the spin den- 
sity in the large dot at the tunneling position r = and 

We also consider (see Sec. ITTJ) a multi-orbital "Anderson"- 
type model that explicitly includes the effect of charge fluctu- 
ations on the quantum dot S: 

-^int = E m4w ( ^ii» + tm[fo a d ma + H.C.] 

ma ma 

+ U(N s -n s g ) 2 . (3) 

Here the quantum dot S is described by a set of spin- 
degenerate energy levels ej^ created by d^ which couple 
to the state /q ct = ^2 a (j) a {Qi)c i a(T in R. Interactions are in- 
cluded through the usual charging term of strength U, where 
N s = J2ma dma-dma and is the dimensionless gate voltage 
applied to the small dot. When m takes only a single value, 
this reduces to the usual single-level Anderson model. The 
crucial feature of this model is that the R-S tunneling term 
(proportional to i) involves only one state in the reservoir. 

For temperature T much larger than not only the mean level 
spacing Ar but also the corresponding Thouless energy of 
the reservoir dot, the discreetness of the spectrum as well as 
mesoscopic fluctuations in R can be ignored. Thus one ex- 
pects to recover the traditional behavior of a spin- 1/2 Kondo 
or Anderson model. If T<g; Ar, however, significantly differ- 
ent behavior is expected. A simplifying feature of this limit 
is that many physical quantities can be derived simply from 
properties of the ground state and low-lying excited states. 

To study the low temperature regime, we shall therefore in 
a first stage consider the low energy (many-body) spectrum 
of the Hamiltonian Eq. (Q}. Specifically, in Sec. [H] we extend 
(slightly) a theorem from Mattis 43 that enables us to infer the 
ground state spin of the system. Using weak and strong cou- 
pling perturbation theory, we then construct in Sec.|III]the fi- 
nite size spectrum of the Kondo problem in a box. 

In a second stage, we consider a few observable quantities 
that are derived simply from the low energy spectra. We start 
in Sec. [IV] with tunneling spectroscopy, obtained by weakly 
connecting two leads to the reservoir dot (Fig.[T|i. Using a 
rate equation approach, we predict generic features in the non- 
linear I-V of our proposed device. We then address in Sec. [V] 
the low temperature magnetic response of the double dot sys- 
tem, and in particular discuss the difference between local and 
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FIG. 2: (Color online) Mapping of the Kondo problem into a spin- 
chain with an impurity. The site '0' is the point in the bath that 
interacts with the impurity, and is used as the site to begin the tri- 
diagonalization of the one-body Hamiltonian of the bath. The sites 
are labeled by the integer i, in the sequence they are generated by the 
tri-diagonalization procedure. After tri-diagonalization, we are left 
with non-interacting electrons that feel an on-site potential oti and 
can hop only to the neighboring sites with amplitude ti^+i. 



impurity susceptibilities which, although essentially identical 
for T > A R , differ drastically when T < A R . A further is- 
sue that we study is that the charging energy in R fixes the 
number of electrons rather than the chemical potential; thus, 
the canonical ensemble must be used rather than the grand- 
canonical. Use of the canonical ensemble accentuates some 
features in the susceptibility. Finally, we conclude in Sec. lVIl 



H. GROUND STATE THEOREM 

We now prove an exact ground state theorem for the mod- 
els defined in Eqs. <[TJ-<[3j : the ground state spin of the system 
is fixed, and in particular cannot depend on the coupling be- 
tween the small dot and the reservoir. We give the value of 
this ground state spin in a variety of cases. 

The theorem is mainly an extension of a theorem due to 
Mattis. 43 It relies on the fact that in a specially chosen many- 
body basis, all the off-diagonal matrix elements of these 
Hamiltonians are non-positive. It is then possible to invoke 
a theorem due to Marshall 44,45 to infer the ground state spin - 
a proof of Marshall's sign theorem is in AppendixlAl 



A. "Kondo"-type models 

The starting point of the proof is the tri-diagonalization of 
the one-body Hamiltonian of the reservoir, Hr. Beginning 
with the state / U(T , one rewrites Hr as a one dimensional chain 
with only nearest-neighbor hopping. 1,46 This transformation is 
illustrated in Fig. [2] In this one-body basis, the "Kondo"-type 
model Eqs. ([Hi-© can be rewritten as a sum of a diagonal and 
off-diagonal part: 



Hot, = -!^[S+«-(0)+H.c.] 



(6) 
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Condition 1 of the the Marshall theorem requires us to find a 
many -body basis in which all the off-diagonal matrix elements 
are non-positive. Consider the following basis: 

= (-ir- 5 /^.../^,.../* |0> ® |m> (7) 



with m the quantum number of S z of the local spin and the 
site labels (positive integers) ordered according to i% < . . . < 
i N R and j\ < . . . < j N n . Note that this basis is diagonal 

with respect to both the total number of electrons in R, N R — 
N R + N^, and the z-component of the total magnetization, 

Sf ot =m+(iVf-iVf)/2. 

The off-diagonal matrix elements come from two terms, 
the spin-flip term and the fermion hopping. With regard to 
the fermion hopping term, first, since the fermions have been 
written as a one-dimensional chain, there is no sign from the 
fermionic commutation relation. Additionally, one can use 
the freedom to choose the phase that defines the one-body 
states fj a to make the hopping integrals i^j+i negative. Since 
the number of phases is the same as the number of hopping 
integrals ti,{+i, all the tj,j+i can be made negative, as in 
Eq. ©. This ensures that all off-diagonal matrix elements of 
the fermion hopping term in the many-body basis, Eq. @, are 
non-positive. With regard to the spin-flip term, note that its 
sign in iJoD can be fixed by rotating the spin S by an angle 
7T about the z-axis. In order to ensure that the off-diagonal 
elements due to the spin-flip term are negative, we have to in- 
clude the additional phase factor (— l) m ~ 5 appearing in the 
definition of the basis states in Eq. (0. 

Since the basis Eq. ((T) is diagonal in Nr and S* ot , we will 
work in a fixed (Nr, S* ot ) sector. Condition (ii) of Marshall's 
theorem - connectivity of the basis states by repeated appli- 
cation of Hjis - is easily seen to be satisfied for the "Kondo" 
model for all J ^ 0, in a given (Nr, S* ot ) sector. However, 
when J = 0, condition 2 is violated: the impurity spin cannot 
flip and hence some basis states in a (Nr, S z ot ) sector cannot 
be connected to each other by repeated applications of Hrs. 

We have thus shown that the Kondo model satisfies the two 
conditions of Marshall's theorem in a given (Nr, S z ot ) sec- 
tor. Now note that given Nr and S, the competing spin mul- 
tiplets for the ground state spin (Stot) can either be integer 
spin multiplets or half-integer spin multiplets. Suppose for in- 
stance they are integer multiplets (this is true, e.g., when Nr 
is odd and S = 1/2). Marshall's theorem guarantees that in 
the S z ot = sector the lowest eigenvalue can never have a 
degeneracy; this ensures that in a parametric evolution there 
can never be a crossing in the S^ ot = sector. Since each 
competing multiplet has a representative state in the S* ot — 
sector, we infer that the ground state spin does not change as 
the coupling J is tuned. This is true as long as we do not 
cross the point J — 0, because this point (as explained above) 
violates condition 2 in the proof of the theorem. Hence the 
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TABLE I: Ground state spins for different Kondo problems according 
to the theorem combined with perturbation theory. Marshall's theo- 
rem adapted to the model defined by Eqs. (Q3 and l(2j says that the 
ground state spin does not change in a parametric evolution of the 
Hamiltonian. The only exception is the crossing of the point J = 0, 
hence the sign of J appears in the table. 

ground state spin can be different for ferromagnetic and anti- 
ferromagnetic J, but does not change with the magnitude of 
the coupling: the ground state spin for all J may hence be in- 
ferred by lowest order perturbation theory in J. The ground 
state spin for a few representative cases is displayed in TableJ] 

B. "Anderson"-type models 

We can prove a similar theorem for the model defined by 
Eqs. (fl~|i and (01. We begin by tri-diagonalizing the electrons 
in the reservoir R, as for the Kondo case. In addition we have 
to tri-diagonalize the electrons in the quantum dot S, a process 
which begins with the state d 0a = (1/£r,-s) J2 m tm<$mo where 
tR-s = (J2 m tm) 1 ^ 2 - An organization of H R „ S into diagonal 
and off-diagonal parts then yields 

HD = "m dLJma + U(N S - n S g ) 2 (8) 

m,<7 

Hqd = -|tR-s|(/ t CT 4r + H.c.) 

_ y^(l*m,m+ll d^dm+la + H.C.) 

fiafi+lcr + H.C.) . (9) 

We note again that the sign of all the hopping integrals can be 
fixed as displayed above by an appropriate selection of the ar- 
bitrary phase that enters the definition of the and the fj a . 
The appropriate basis that has only non-positive off-diagonal 
matrix elements is, then, simply 

T T 

4i-^i4i-4aJ°> ■ (10) 
l I 

The total number of particles is now N tot = Nf + Nf + + 
Nf, and the z-component of spin is Sf ot = (Nf +Nf- Nf- - 
Nf)/2. 
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FIG. 3: Weak-coupling perturbation theory: schematic illustration of 
the unperturbed system for (a) N R odd and (b) N R even. The spin 
marked S is that of the small quantum dot while the solid lines rep- 
resent the spectrum of the (finite) reservoir. The dashed lines show 
the lowest energy orbital excitation in each case; note that in the even 
case, any excitation requires promoting an electron to the next level 
and so involves a minimum energy of order Ar. 

In the case of the Anderson-type model Eqs. (fTJ and (0, 
the result for the ground state spin is remarkably simple: the 
ground state spin has 5 tot = for iV tot even and 5 tot = 1/2 
for 7V tot odd. There is no possibility of having a ground state 
spin other than the lowest. 

III. FINITE SIZE SPECTRUM 

In this section, we outline the main features of the low- 
energy finite-size spectrum for the Kondo problem, Eqs. (fl~|i- 
(O. The basic idea is to use perturbation theory around its 
two fixed points: at the weak coupling fixed point ( J = 0) ex- 
pand in J, and at the strong coupling fixed point expand in the 
leading irrelevant operators (Nozieres' Fermi-liquid theory). 
We begin by analyzing the classic case of S = 1/2 with anti- 
ferromagnetic coupling, and then turn to the under-screened 
Kondo problem realized by anti-ferromagnetic coupling and 
5 = 1. 



A. 5 = 1/2: Screened Kondo problem 

Weak-Coupling Regime: In the weak coupling regime de- 
fined by Ar 3> Tk, given a realization of the reservoir R, we 
can always make J small enough so that the spectrum can be 
constructed through lowest order perturbation theory. 

The unperturbed system for N odd is shown schematically 
in Fig.[3}a). At weak coupling the eigenstates follow from 
using degenerate perturbation theory in all the multiplets of 
the unperturbed system. The ground state and the first excited 
state are obtained by considering the coupling 

#top = </|0a top (O)| 2 S top -5 (11) 

where s top is the spin of the topmost (singly occupied) level 
atop of the large dot. The ground state is therefore a singlet 
( J > 0) and the first excited state is a triplet with excitation 
energy 

SE ST = J|^ top (0)| 2 «A R . (12) 



FIG. 4: Strong-coupling perturbation theory: schematic illustration 
of the unperturbed system for (a) 7V R odd and (b) N R even. For 
5 = 1 /2, a conduction electron is bound to the impurity at strong 
coupling. This leaves effectively a gas of even (odd) weakly interact- 
ing quasi-particles. We note here that in the usual case of Tk <S D 
(bandwidth of the reservoir), the formation of a singlet between the 
impurity spin and the reservoir for Tk 3> Ar is a complicated many- 
body effect: it is no? just a singlet between the spin and the topmost 
singly occupied level, as is the case for very weak coupling. 

The next excited states are obtained by creating an electron- 
hole excitation in the reservoir [shown as a dashed arrow in 
Fig. [3fa)] . Combining the spin 1 /2 of the reservoir with that of 
the small dot, one obtains a singlet of energy ~ Ar separated 
from a triplet by a splitting ~ J'Ar /4, where we define J = 
J p with p — (|0 q (O)| 2 )/Ar the mean local density of states. 

In the N even case depicted in Fig.|3b), the ground state 
is trivially a doublet. The first excited eigenstate of the un- 
perturbed system is an 8-fold degenerate multiplet obtained 
by promoting one of the bath electrons to the lowest available 
empty state [shown as a dashed arrow in Fig. [2b)]. As J is 
turned on, this multiplet gets split into two Stot = 1/2 doublets 
and one Stot = 3/2 quadruplet. In general the two doublets 
have lower (though unequal) energy than the quadruplet. 

Strong-Coupling Regime: For Tk S> Ar , on the other hand, 
the impurity spin S is screened by the conduction electrons, 
and we can use Nozieres' "Fermi-liquid" theory. 47 ' 48 In the 
very strong coupling limit, one electron is pulled out of the 
Fermi sea to bind with the impurity; this picture essentially 
holds throughout the strong coupling-regime. 4,47,48 For N odd 
(even) one ends up effectively with an even (odd) number of 
quasi-particles that interact with each other only at the impu- 
rity site through a repulsive effective interaction 

f/ FL ~ (A2/T K )n T (0)n ; (0) (13) 

which is weak (Tk 3> Ar). The quasi-particles have the same 
mean level spacing Ar as the original electrons, but the spac- 
ing between two quasi-particle levels is not simply related to 
the spacing of the original levels in the chaotic quantum dot. 
This case is illustrated in Fig.|4] 

For TV odd, the ground state is thus a singlet (as expected 
from our theorem), and the excitations start at energy ~ Ar 
since a quasi-particle must be excited in the reservoir. The 
first two excitations consist of a spin S to t = 1 an d a Stot = 
0. Because the residual quasi-particle interaction is repulsive, 
the orbital antisymmetry of the triplet state produces a lower 
energy; the splitting is about ~ A^/T K . 

In the N even case at strong coupling, there are an odd num- 
ber of quasi-particles in the reservoir, and so the ground state 



J J 

a) N = odd b) N = even 

FIG. 5: (Color online) Schematic of the energy eigenvalues of the 
double-dot system as a function of the coupling J for (a) N odd and 
(b) TV even in the 5 = 1/2 anti-ferromagnetic coupling case. En- 
ergy differences are shown with respect to the ground state which 
therefore appears on the x-axes. The relation to the double-dot ex- 
periment Fig.[T]is that the y-axis here is like Vbias and the x-axis is 
like Vpinch- The excitations will show up as peaks in the differential 
conductance G. 

is a doublet. The first excited multiplet must involve a quasi- 
particle-hole excitation in the reservoir. There are two such 
excitations that involve promotion by one mean level spacing 
on average (either promoting the electron in the top level up 
one, or promoting an electron in the second level to the top 
level). Thus, the first two excitations are doublets. 

Crossover between Weak- and Strong-Coupling: Remark- 
ably, the ordering of the Stot quantum numbers of the ground 
state and two lowest excitations is the same in both the T K 3> 
Ar and T K <C Ar limits. It is therefore natural to assume that 
the order and quantum numbers are independent of T k /Ar. 
Thus we arrive at the schematic illustration in Fig. [5] 

B. 5 = 1: Under-screened Kondo problem 

The theorem and perturbation theory analysis presented 
above has an interesting generalization to the under-screened 
Kondo effect, in which S > 1/2. We will consider for con- 
creteness the case 5 = 1. Note that the under-screened Kondo 
effect has been realized experimentally in quantum dots42 

For N even and S = 1, we find that the ground state for 
all J has Stot — 1- At weak coupling, this follows directly 
from perturbation theory - the reservoir has spin zero and J 
is too small to promote an electron so the spin of the ground 
state is just that of the small dot. The theorem then implies 
that S = 1 for all J. The first excited multiplet is at energy of 
order Ar. It splits into a singlet, two triplets, and a quintuplet; 
as J increases, the singlet has the lowest energy because the 
coupling is anti-ferromagnetic. 

In the opposite limit of strong coupling, as J — * oo, one 
electron from R binds to the impurity spin forming a spin- 1/2 
object. For J = oo, this spin does not interact with the quasi- 
particles in R; however, when J ^ oo, the flow to strong cou- 
pling generates other irrelevant operators that connect the spin 
to the quasi-particles. It is known from studies of the under- 
screened Kondo problem that the leading irrelevant operator 
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a) N = even b) N = odd 

FIG. 6: (Color online) Schematic illustration of the finite-size spec- 
trum of the under-screened Kondo problem, 5 = 1 with anti- 
ferromagnetic coupling, for (a) TV even and (b) TV odd. In each panel 
we show all the excitations up to order A, both at strong and weak 
coupling. When it seems plausible, we have connected the strong 
and weak coupling limits; note the necessity of crossings among the 
excited states in (a). 



is a ferromagnetic Kondo coupling (the sign of the coupling 
follows heuristically from perturbation theory in t/J). How- 
ever, since one of the electrons is bound to the spin, there is an 
odd number of quasi-particles in the effective low energy fer- 
romagnetic Kondo description - the level filling is as in Fig. 
Htb). Since the ferromagnetic Kondo problem flows naturally 
to weak couplings ' we are again justified in doing perturba- 
tion theory in the coupling, and so recover that the ground 
state has Stot = !■ From the small ferromagnetic coupling, 
we conclude that the first excited state is a singlet separated 
from the ground state by an asymptotically small energy (as 
Tk/Ar — > oo). The next excited state involves promotion of 
a quasi-particle to the next level within R and so has energy 
of order Ar. It is a triplet because of the ferromagnetic cou- 
pling, with a nearby singlet in the strong coupling limit. Note 
that there are two possible quasi-particle excitations with en- 
ergy of order Ar (as discussed in the 5 = 1/2 case), and so 
two singlet-triplet pairs. 

The proposed crossover from weak to strong coupling for 
TV even is shown in Fig.[6ja). Note that in this case, level 
crossings of excited states must occur: the two singlets at 
energy of order Ar at strong coupling come from energies 
greater than Ar at weak coupling, and so cross the S= 2 state. 
The two singlet-triplet pairs at strong coupling are shown to be 
slightly different because each involves a different level spac- 
ing; thus, there is an additional level crossing as one of the 
singlets comes below a triplet. 

In the TV odd case, weak anti-ferromagnetic coupling im- 
plies that the ground state spin is Stot = 1/2. The first ex- 
cited state is the other multiplet involving no excitations in 
the reservoir, Stot = 3/2. The next excited states are the 
Stot = 1/2 and 3/2 states that involve promoting one elec- 
tron by one level. In the strong coupling limit, we repeat 
the mapping to a ferromagnetically coupled impurity, yielding 
this time an even number of quasi-particles in the reservoir. 
Now the first excited state involves promoting a quasi-particle 



in the reservoir by one level; the ferromagnetic coupling im- 
plies that the Stot = 3/2 state has the lowest energy among the 
possible multiplets. Making again the reasonable assumption 
that the two limits are connected to each other in the sim- 
plest manner possible, we arrive at the schematic illustration 
in Fig.|6jb). In contrast to the TV even case, no level crossings 
are definitely required. 

We stress here that the evolution of the finite size spectra 
shown in Figs.[5]and[6]are totally different in each of the cases 
illustrated. The finite-size spectrum is hence an interesting 
way to observe the Kondo effect in nano-systems, each impu- 
rity problem having its own unique spectrum. 



IV. NON-LINEAR I V CHARACTERISTICS OF THE R-S 
SYSTEM 

We now turn to the question of how to observe the features 
of the finite size spectrum delineated in the previous section. 
Any physical observable depends, of course, on the spectrum 
of the system and so could be used as a probe. We choose 
to concentrate on two: (1) In the next section, we discuss the 
magnetic susceptibility of the R-S system, a classic quantity in 
Kondo physics. (2) In this section we discuss the conductance 
across the device shown in Fig.Q] The advantage of this phys- 
ical quantity is that the finite-size spectrum can be observed 
directly in the proposed experiment. The emphasis here is on 
transfer of electrons entirely by real transitions; cotunneling 
processes, which involve virtual states, are briefly discussed 
at the end of the section. 

A current through the R-S system clearly involves number 
fluctuations on it. For a general value of the gate voltage [n^ 
in Eq. ([T)], however, the ground state will have a fixed num- 
ber of electrons, and hence G = (Coulomb blockade). When 
Vbias = V1 — V2 is increased sufficiently, the Coulomb block- 
ade is lifted, and G(Vbias) nas a sequence of peaks. It is 
possible to extract the excitation spectrum of the R-S system 
from the position of these peaks. 51 In principle, there is a peak 
in G for every transition a — > (3 that involves a change in 
TV. As discussed in subsection llV Cl we shall however choose 
a particular setting such that only a limited number of these 
transitions play a role, making in this way simpler the recon- 
struction of the underlying low-energy many-body spectra. 



A. Method 

In order to describe transport through the R-S system (re- 
alized through either a double dot or a metallic grain with a 
single magnetic impurity), we solve the appropriate rate equa- 
tions for the real transitionsi 5 -^ The rate equations are a limit 
of the quantum master equation in which the off-diagonal el- 
ements of the density matrix are neglected. The dynamics of 
the quantum dot can then be described simply by the proba- 
bility P a that the R-S system is in a given many-body state a. 
In thermal equilibrium these P a are the Boltzmann weights. 
The electrons in the leads are assumed to always be in ther- 
mal equilibrium; hence, the probability that a given one-body 
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state in the leads is occupied is given simply by the Fermi- 
Dirac function /(e) = l/(e e / T + l). Here, e is the deviation 
from the electro-chemical potential Ep + V12, where V\ and 
V2 are the voltages on leads 1 and 2 respectively. 

Steady state requires that the P a are independent of time. 
Hence, the various rates of transition from a to /?, Kp a , must 
balance, leading to a linear system for the P a , 

2^ = 2^. (14) 

In addition, the occupation probabilities should be normal- 
ized, J2a P a = l - 

There are four transitions that have to be taken into account: 
addition or removal of an electron from lead LI or L2. We 
denote the rates for these four processes as A^ 1,L2 , and the 
A a j3 in Eq. ([Pil l are sums of these four transition rates. Once 
we have the P a from ( TT4l i, the current is simply 

= = D< - a ^ L2)Pq ■ (15) 

The conductance G then follows by differentiating /(Vbias)- 
The rates A a p can be calculated in second order perturba- 
tion theory in the reservoir-lead coupling term, using Fermi's 
golden rule. 51 For example, consider the addition of an elec- 
tron to R-S from Li corresponding to a transition /3(N) — > 
a(N + 1) on R-S: 

= Jdep(e)f(e)S(E a -E fj -e-V 1 ) 
= F 1 f(E a -E -V 1 ) (16) 

where Ti^ = 2irVf 2 p(Bf) / 'h. V\ is the amplitude for the 
above process. Although, in general it will have some depen- 
dence on a and j3 as well as the coupling J, we will ignore 
such dependence here. We will, however, retain the distinc- 
tion between Vi,2 and allow these to be tuned by the gates 
that define the R-Ll and R-L2 junctions. 

To summarize our approach, to find the conductance in the 
proposed tunneling experiment, we have solved the rate equa- 
tions for transferring an electron from lead 1 to the reservoir 
and then to lead 2^!^ We assume that (1) the coupling of the 
lead to each state in R is the same (mesoscopic fluctuations 
are neglected), (2) the Kondo correlations that develop in R-S 
do not affect the matrix element for coupling to the leads, 53 
(3) there is a transition rate A re i that provides direct thermal 
relaxation between the eigenstates of R-S with fixed N, (4) 
the electrons in the lead are in thermal equilibrium, and (5) 
the temperature T is larger than the widths Y\, T2 of the R-S 
eigenstates due to LI and L2. 



B. Magnetic Field 

A Zeeman magnetic field Bz can be used as an effective 
probe of the various degeneracies of the R-S system. We shall 



assume that the magnetic field does not couple to the orbital 
motion of the electrons: 

H z = -SMbBz • S . (17) 

This can be achieved in the semiconductor systems by apply- 
ing the field parallel to the plane of motion of the electrons. 
The effect of an orbital magnetic field in ultra-small metallic 
grains is argued to be small in Ref. JHfor moderate fields. 

We may neglect the effect of Bz on the lead electrons: The 
only characteristic of the lead electrons appearing in the rate 
equation calculations is the density of states at Ep. All that Bz 
does to the lead electrons is to make the modification p(Ep) — > 
p(Ep ± gfisBz/2). Since the band is flat and wide (on the 
scale of Bz) to an excellent approximation, this has no effect. 

The effect of Bz on the R-S system is complicated if the g- 
factors for the S and R electrons are different, as would be the 
case for a magnetic impurity in a metallic nano-particle. If we 
assume, however, that the g-factors for the electrons on the S 
and R systems are the same, as is relevant for the semiconduc- 
tor quantum dot case illustrated in Fig.Q] then Hz becomes 
simply —gfXBB%S? ot . The energy of a given many -body level 
a is then E a — gp,BBzS^ ot where S* ot is the corresponding 
eigenvalue of the many-body state. 

C. Application to R-S System 

To identify characteristic features in the transport proper- 
ties, let us analyze a situation in which only a limited number 
of transitions show up. 51 For a S = 1/2 Kondo problem the 
most interesting features appear in the spectrum when there is 
an odd number of electrons in the reservoir. These states ap- 
pear clearly when an electron is added to a A even reservoir 
and the parameters are such that the excited states of the N+l 
electron reservoir dominate. 

We thus consider the following situation: For zero bias, as- 
sume that the R-S system is brought into a Coulomb blockade 
valley, not far from the N — ► N + 1 transition. This could 
be done by adjusting V\ and V2 in the setup of Fig.Q] (with 
Vi = V2), or more realistically with the the help of the addi- 
tional gate voltage in Eq. (fl}. We take this setup as the 
origin of the bias potentials (Vi = V% = 0). Upon applying a 
bias V\ > V2, electrons flow from lead 1 to lead 2. 

We assume that the rates T12 are sufficiently small that vir- 
tual processes (cotunneling) can be entirely neglected for this 
subsection; that is, all relevant transitions occur on shell and 
can be described by the Fermi golden rule expression Eq. ( TTol ). 
Furthermore, we take T2 Ti, Because Vi > V2, this means 
that it takes much longer to add an electron to the dot than to 
empty it. Thus, the dot tends to be occupied by N electrons. 

Several conditions are needed in order to restrict the dis- 
cussion to just the lowest lying states of the system. First, 
we shall assume that T<Aj; so that in equilibrium only the 
ground state 5=1/2 doublet, with energy Eg [N], needs to be 
considered. In a non-equilibrium situation, however, higher 
excited states Ei[N] ~ A# can also be populated: the ex- 
cess energy of the electron supplied by the bias can be used 
to leave the dot in an excited state. Apart from the ground 
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BIAS 

FIG. 7: (Color online) Illustrative calculation of transport spec- 
troscopy starting from the ground state of a N even Kondo R-S sys- 
tem. Transitions are S = 1/2 — > S = 0, 1 as marked. Top panel: 
The differential conductance as a function of bias voltage for differ- 
ent values of the asymmetry between L, R tunneling rates. Bottom 
panel: The probability of occupation of the two states forming the 
S = 1/2 ground state. For large asymmetries P| + Pj. fa 1, as ex- 
pected. Even though ksT / gusBz =0.2 <C 1, P± is large because 
this calculation neglects inelastic relaxation on the R-S system, al- 
lowing it to stay well out of equilibrium. T\ = 0.01, gfisBz = 0.4, 
and T 2 /Ti = 10, 3, and 1 from top to bottom.— 




BIAS 

FIG. 8: (Color online) Transport spectroscopy with energy relaxation 
included. I^/Fi = 10 (fixed), the parameter A ro i that models energy 
relaxation varies, and other parameters are as in Fig.|7] Top panel: 
The differential conductance as a function of bias voltage. Note how 
the symmetry of the 5=1/2 — > S = 1 peaks is unaffected while 
the S — 1/2, S z = —1/2 S = transition is suppressed (this 
transition would completely vanish if the R-S system were in thermal 
equilibrium). Bottom panel: The probability of occupation of the two 
states forming the 5=1/2 ground state. As A rc i is increased, the R-S 
system has a larger probability to occupy its ground state (S z = 1/2). 



state doublet, we take all excited states Ei [N] to be higher in 
energy than E G [N + l]£z Then, if an excited state is pop- 
ulated, it will quickly relax to an energy below Eq[N + 1] 
through a rapid exchange of particles back and forth between 
the dot and lead 2 (1^ 3> Ti). Because off-shell processes are 
assumed negligible, this relaxation will stop as soon as an N- 
electron state Ei [N] below Eg [N + 1] is reached. We assume 
that the energy of the first excited 5=1/2 doublet, E\ [N], is 
large enough that E X [N\ + V 2 > E G [N + 1} > E G \N] + V%. 
Then only the TV-electron ground state multiplet needs to be 
retained in the calculation. 

With regard to the N+l electron states, we limit ourselves 
to a small enough bias such that only transitions to the three 
lowest excited multiplets need to be taken into account. In this 
way, only a small number of transitions will show up in the 
excitation spectrum, making it relatively simple to analyze. 51 

When a magnetic field is applied, note the following un- 
usual behavior: Since there is no way to decay from the S z = 
- 1/2 state to the S z = +1/2 state of the lowest doublet with- 
out involving virtual processes explicitly neglected here, the 
doublet will remain out of equilibrium: the S z = — 1 /2 state 
can be significantly populated even though g/isBz ^S> ksT. 



D. Results for G 

With the above assumptions, results for the differential con- 
ductance are shown in Fig. [7] We assumed that the system pa- 
rameters (gate potential, Vi, V2, and SEst) are such that the 



first three states of the N+l electron system coincide with the 
ground state of the N electron system for Vbias = 1> 2, and 3, 
respectively, at B = 0. (See Fig 4. of Ref.Hlfor dl/dV in the 
B = case.) We are thus assuming that the excited triplet state 
lies midway between the two lowest singlet states (see Fig.|5]l, 
placing ourselves in the middle of the cross-over regime. 

First, note that the ground state to ground state transition, 
5=1/2 — * 5 = 0, yields only one peak even at non-zero 
Bz- This is because the S z = — 1/2 state of the doublet cannot 
be populated before some current is flowing through the R- 
S system {k-&T / g^sBz = 0.2 <§C 1). However, after the first 
transition, the N+l state can decay into the S z — — 1/2 state. 
Hence, we expect the higher 5=1/2 — > 5 = transitions to 
split in a magnetic field, as in Fig. [7] 

The next feature to understand is the two 5=1/2 — * 5 = 
1 transitions. These two peaks occur because out of the six 
transitions between the multiplets, two are forbidden by spin 
conservation and the other four split into two degenerate sets. 

How is one then to distinguish between a 5 = and 5=1 
state, since they both split into two as a function of B%! One 
possible method is to observe the peak heights in G( Vbias) 
keeping Bz fixed. These are plotted in Fig.|7]for a variety of 
r 2 /ri^ A clear feature is that the two 5 = 1/2 — > 5 = 
peaks are very asymmetric, while the 5 = 1/2 — * 5 = 1 are 
almost symmetric. This is for a robust physical reason: each 
5=1 peak gets contributions from both S z = l/2 and —1/2 
initial states, while in the 5 = transition each peak gets a 
contribution from only one, the S z = 1/2 for the taller peak 
and S z = —1/2 for the shorter one. The associated proba- 
bilities, Pf and Pj, are shown in the lower panel of Fig. [7] 
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Thus the peak heights in the 5 — 1 transitions are insensitive 
to the difference between the probability of occupation of the 
two states in the doublet, while the peak heights in the 5 — 
transitions are sensitive to this difference. 



E. Energy Relaxation 

In any real system, there are mechanisms of energy re- 
laxation beyond the energy conserving exchange of electrons 
with the leads that is given by the rate equations. These mech- 
anisms can involve, for instance, interactions with phonons or, 
more simply, higher order virtual processes between the R-S 
system and the leads that are neglected in the Fermi's golden 
rule approach Eq. <[T6t . These relaxation processes are partic- 
ularly important for the second 5 = 1/2 — > 5 = transition. 
If the system is in perfect thermal equilibrium this transition 
should yield a single peak, even in the presence of a Bz ^ 0. 
The second peak is suppressed even if only on-shell processes 
are taken into account, as discussed above, but explicit energy 
relaxation causes this suppression to be more pronounced. 

To model energy relaxation, we include a transition rate be- 
tween the S z = ±1/2 states that satisfies detailed balance (i.e., 
with Boltzmann weights), 



,s a /T 



-e a /T + e -ep/T 



(18) 



where e a is the energy eigenvalue of the a ih state. 

The effect of this term is shown in Fig. [8] Clearly as 
the relaxation rate is increased, the peak in the second 5 = 
1/2 — > 5 = transition coming from non-equilibrium effects 
is suppressed further. Note, however, that the heights of the 
5=1/2 — * 5=1 transition are unaffected (in both relative 
and absolute magnitude). 



F. Cotunneling Spectroscopy 

While the approach proposed above should be reasonably 
simple to implement, because the excitations of both the N 
and N + 1 electron systems may come into play, the resulting 
experimental conductance curves may in some circumstances 
be non-trivial to interpret. Therefore, we mention, without 
going into detail, an alternative way to extract the excita- 
tion spectra from the differential conductance. Though within 
the "Coulomb blockade diamond", on-shell processes such as 
the ones considered above are forbidden by energy conserva- 
tion constraints, a small current can nevertheless be measured, 
which is associated with virtual (cotunneling) processes.— 

At very low bias, these virtual processes are necessarily 
elastic as the electron transferred from one lead to the other 
does not have enough energy to leave the R-S system in an 
excited state. However, each time Vbias reaches a value 
corresponding to an excitation energy of the system with N 
electrons, a new "inelastic" channel is open, as the electron 
has the option to leave the R-S system in an excited state 
as it leaves the structure. The opening of these new chan- 
nels produce steps in the differential conductance within the 



Coulomb diamond. These steps are small, but clearly observ- 
able experimentally i 57 ' 58 ' 59 

Because of the smallness of the associated currents, observ- 
ing this substructure within the Coulomb diamond is certainly 
more challenging experimentally than for observing the main 
peaks associated with on-shell processes. On the other hand, 
the time elapsed between the successive transfers of an elec- 
tron across the structure is large enough that the initial state 
of the N -particle system is always the ground state. If they 
can be measured accurately, the cotunneling steps within the 
Coulomb diamond may therefore lead more directly to the TV- 
particle excitation spectra. 

Summarizing this section, we have shown in detail how 
dl/dV measurements enable one to extract the finite size 
spectrum and spin quantum numbers of the R-S system, us- 
ing the case when the ground state has an even number of 
electrons as an example. In particular, we argued that the rel- 
ative peak height of the Zeeman split terms {Bz ^ 0) reflects 
the spin quantum number of the excitation: asymmetric peak 
heights correspond to 5 = 0, whereas symmetric peak heights 
correspond to 5 = 1. The case when N is odd is straightfor- 
ward to analyze in a similar way. Transitions from 5 = to 
5= 1/2 or 3/2 can easily be distinguished: the former splits 
into two in a magnetic field while the latter splits into four. 



V. MAGNETIC RESPONSE OF THE 
DOUBLE DOT SYSTEM 

We turn now to studying a second physical observable 
which probes the finite size spectrum of the system, namely 
the magnetic susceptibility defined by 



X = s 



1 d 2 \ogZ 
dB 2 



(19) 



where Z is the canonical or grand-canonical partition func- 
tion depending on the ensemble considered. As in Section 
IIVI we assume that the magnetic field is in plane so that only 
the Zeeman coupling needs to be considered, Eq. (fTTT i. We 
furthermore distinguish between the local susceptibility Xioc, 
corresponding to the case where B couples only to the quan- 
tum impurity spin (S = S), and the situation where B couples 
to the total spin of the R-S system (S = Stot)- In the latter 
case, the impurity susceptibility Ximp is defined as the differ- 
ence Xtot - Xo between the total magnetic response and that of 
R in the absence of the impurity dot. 

For a wide (D > T K ) and flat (T > E* h ,A R ) band the 
local and impurity susceptibilities are essentially identical. 60 
Indeed the effect of the magnetic field on the reservoir elec- 
trons is just to shift the energies of the spin up electrons with 
respect to the spin down by a fixed amount. If the spectrum is 
featureless, this only affects in practice the edge of the band, 
which in the limit J = J p — > and D — > oo with fixed Xk 
will not affect the Kondo physics. More precisely, for small 
but finite J p (and again for a wide flat band) the impurity sus- 
ceptibility, being associated with the correlator of a constant 
of the system, can be written as 



f x (T/T K ) + 0(T K /D) 



(20) 
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where f x is a universal function of the ratio (T/Tk). On 
the other hand, since the spin S of the impurity is not con- 
served, a multiplicative renormalization factor z Xloc needs to 
be introduced for the local susceptibility so that Tkxioc — 
z Xloc / x (T/Xk). For z Xlac we use a form motivated by two 
loop renormalization, l/z Xloc = 1 — J + a J 2 for J -C 1, 
with the coefficient of the quadratic term determined empiri- 
cally, a — —0.4. (For a discussion of z Xloc in the context of 
two loop renormalization, see e.g. Ref.|6j].) We note here that 
in the universal regime T^/D — > 0, one also has J — > and 
hence z Xloc — > 1. In practical numerics, even though T^/D is 
small enough that the 0(Tr/D) correction can be neglected, 

— l/ln(D/Tk) need not be as small; hence, it is neces- 
sary to include the prefactor correction z Xloc to observe good 
scaling behavior. 

In the regime T <C Ar that we consider here, however, the 
reservoir electron spectrum is not featureless near the Fermi 
energy, and the Zeeman splitting of the conduction electrons 
affects in practice the whole band, and not just the band edge. 
One therefore does not particularly expect any simple relation 
between the local and impurity susceptibilities. We now dis- 
cuss the behavior of these quantities in this regime. We start 
with the canonical ensemble, for which the number of parti- 
cles N, and therefore the parity of N, is fixed. We'll consider 
in a second stage the grand canonical ensemble and so neglect 
charging effects in the reservoir [Ec — in Eq. ([]}]; in this 
case the spin degeneracy induces finite fluctuations of the par- 
ticle number even in the zero temperature limit. 



A. Canonical ensemble 

Since S* ot is a good quantum number, the impurity suscep- 
tibility in the canonical ensemble follows immediately from 
the information contained in Fig. [5] i.e. from the knowledge 
of the total spin and excitation energy of the first few many- 
body states. Neglecting all the levels with an excitation energy 
of order Ar (because T<Ar), we simply get for xtot a spin 
1/2 Curie law for even N, and a spin 1 Curie law damped by 
exp[— (3 6 Est] for odd N. In this latter case, the magnetic re- 
sponse in the absence of the impurity is also a spin 1/2 Curie 
law; thus, for 5Est <C Ar one finds 



N even 



Ximp \ 4 

f3{g^ B ) 2 [2eM-P$E S T) - |] , TV odd 

(21) 

The local susceptibility on the other hand involves S which 
is not a conserved quantity. Its computation therefore requires 
knowledge of the eigenstates, in addition to the eigenener- 
gies and total spin quantum numbers contained in Fig. [5] We 
can follow the same approach used in Section|III]and analyze 
the two limiting regimes of coupling between the reservoir 
and the impurity quantum dot. We will then use a numerical 
Monte Carlo calculation in the intermediate regime and in- 
vestigate how well it is described by a smooth interpolation 
between the two limiting regimes. 

In the weak coupling regime, Tr <C Ar, we assume that 



even if some renormalization of the coupling constant J takes 
place, the eigenstates are the ones obtained from first-order 
perturbation theory in this parameter. For even N at T <C Ar, 
the impurity spin decouples from the (frozen) electron sea, 
and one obtains again a spin 1/2 Curie law. For odd N at T<C 
Ar, the system formed by the impurity spin and the singly 
occupied orbital decouples from the set of doubly occupied 
levels. The magnetic response is the same as for two spin 1/2 
particles interacting through Eq. ( fTTb . We thus obtain 



Xlo 
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N even 

iVodd 
(22) 

valid for T K < A R . 

Turning now to the strong coupling regime, we follow 
Nozieres' Fermi liquid picture, 47 ' 48 where low energy states 
(with Ey <C Ik) are constructed from quasiparticles 
which interact locally according to Eq. (JT3J. In a local mag- 
netic field, the energy of a state |W) is modified to 



E*(B) = (gn B B z )(^\S z \^) - (g/isBz) 



2 \ ^ 



(23) 

where the sum is over all the many-body excited states £. The 
first term in this expression yields the effect of a change of the 
quasiparticle phase shift on the energy. It is important when 
N is even: one of the quasiparticle states is singly occupied, 
and its energy is shifted by an amount ~ (g^BBz/Tx) A#. 
Thus the system acts like a spin- 1/2 particle with an effective 
g-factor given by gAn/Tx- The result is a weak Curie sus- 
ceptibility ~ (gfiB^-R/Tfc) 2 /AT at low temperature. 

The second term captures the effect of electron-hole quasi- 
particle excitations. It produces a non-zero contribution even 
when the discreteness of the spectrum is ignored, as may be 
seen as follows. First, the density of states of particle-hole 
excitations of energy AE in a Fermi liquid is proportional to 
AE. In a Kondo state, the density of single-particle states 
is increased by a factor of 1/Tk because of the Kondo res- 
onance. Thus we may replace the sum in Eq. (|23l by an in- 
tegral using a density of states proportional to ~ AE/T K . 
The integral should be cutoff at an energy of order Tk, where 
the Kondo resonance ends. Thus the second term in Eq. ( |23l 
gives a contribution ~ (g/isBz) 2 /Tk to the energy, and a 
corresponding contribution ~ (gfis) 2 /Tk to xioc- 

Note that this second contribution is independent of the 
finite-system parameter A# and so is the universal (bulk) 
part of the local susceptibility. 62 It should behave smoothly 
as T/Ar becomes smaller than one. In particular, if one 
considers a system without mesoscopic fluctuations (i.e. with 
constant spacings and wave function amplitudes at the im- 
purity), we expect to recover the bulk behavior Tkxioc = 
z Xlo J x (T/T K ) for A R /T K -> 0. For N even this relation 
holds only if the weak Curie behavior of the first term is not 
too large; more precisely, we expect \\oc to follow the univer- 
sal behavior as long as T > A 2 r /Tk- 

With these arguments, we have thus arrived at a complete 
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description of the magnetic susceptibility in both the weak and 
strong coupling limits for the canonical ensemble. Note in 
particular the difference in the conditions for having xioc — ► 
Xbuik from those for having Ximp — * Xbulk- Both limits hold in 
the regime T> Ar no matter what the value of Tk- However 
in addition, xioc— >Xbuik for any T as long as Ar/Tr- — > 0. 



B. Grand-canonical ensemble 



Use of the grand canonical ensemble [which involves ne- 
glecting charging effects in the reservoir, Ec = in Eq. ((TJ] 
introduces additional complications compared to the canoni- 
cal ensemble case above. To illustrate, recall first the behavior 
in the absence of the impurity, i.e. for a system of independent 
particles occupying doubly degenerate states e a . For T^> Ar 
in the grand canonical ensemble, the magnitude of the fluc- 
tuation of the number of particles will be significantly larger 
than one. Thus, even if the canonical ensemble result for N 
even is quite different from that for N odd [Eqs. (|2Tt-(l22ll. 
such an odd-even effect would be completely washed out here: 
whatever the choice of the chemical potential \i, configura- 
tions with odd or even N would be as probable. 

In the low temperature regime on the other hand, as soon 
as T <C min Q (|/i — e a \) (which is usually ~ Ar), there is 
a fixed, even number of particles in the system. It is possi- 
ble to make the average number of particles odd by choosing 
[i — e ap for some orbital ap , In that case, as T/Ar — > 0, all 
orbitals a < ap are doubly occupied, all orbitals a > ctF are 
empty, and independent of T the orbital a p has probability 
1/4 to be empty, 1/4 to be doubly occupied, and 1/2 to be 
singly occupied. For quantities showing some odd-even ef- 
fect in the canonical ensemble but no strong dependence on 
TV once the parity is fixed (such as the local susceptibility), 
the grand canonical ensemble produces a behavior which is 
the average of the the odd and even canonical response, even 
though the mean number of particles (N) is odd. 

Turning now to the full R-S system, the above non- 
interacting picture should certainly still hold in the weak cou- 
pling regime. If, either by adjusting fi or by making use of 
some symmetry of the one particle spectrum, (N) is kept fixed 
with an even integer value as T <C Ar, one should recover 
the canonical magnetic response for even N. In contrast, the 
magnetic response for odd (N) should be the average of the 
canonical odd and even responses. 

In the strong coupling regime (following again Nozieres' 
Fermi liquid description), one also has an essentially non- 
interacting picture, but with effectively one less particle since 
one reservoir electron is used to form the Kondo singlet. The 
role of "odd" and "even" are then exactly reversed from the 
weak coupling case: for T <C Ar the grand canonical re- 
sponse will be the average of the canonical odd and even re- 
sponse for even (N), and will be exactly the canonical re- 
sponse for odd (N). 



C. Universality in a clean box 

The previous discussion mainly addressed the two limit- 
ing behaviors - weak and strong coupling. To investigate 
the intermediate regime we now turn to numerical calcula- 
tions. In particular, we use the efficient continuous-time quan- 
tum Monte Carlo algorithm introduced in Ref.|26|, with in 
addition adaptations to compute quantities in the canonical 
ensemble. 63 We study the behavior of the singlet-triplet gap 
5 Est and the local susceptibility xioci the impurity suscepti- 
bility Ximp follows directly from SEst- 

To focus on the consequences of the discreetness of the 
one particle spectrum while avoiding having to explore an ex- 
cessively large parameter space, we disregard the mesoscopic 
fluctuations of the spectrum and the wave-functions. That is, 
we consider the simplified "clean Kondo box" mode l 31 ' 32 ' 33 
defined by cj) a (0) = /oAr and e a+ i — e a = Ar independent of 
a. For initial results for the more realistic "mesoscopic Kondo 
model" see Refs.l25l andl26i 

Under these conditions, the problem is described by only 
three dimensionless parameters: the coupling J = Jp, and the 
two energy ratios D/Ar and T/Ar. For small J and large 
D (3> Ar, T), J and D can be scaled away in the usual man- 
ner so that, except for renormalization prefactors such as z Xlac 
which may still contain some explicit dependence on J, phys- 
ical quantities depend on J and D only through the Kondo 
temperature Tr. 

We therefore expect, again up to the factor z Xloc , that both 
susceptibilities for the "clean Kondo box" model will be uni- 
versal functions of the two parameters T/Ar and Tr/Ar. 
This function may, however, be different for the local and im- 
purity susceptibilities, and will also depend on the parity and 
type of ensemble considered. 

Before discussing how well the limiting behaviors dis- 
cussed above describe the whole parameter range, we shall 
first check that our numerics confirm the expected universal- 
ity. For the impurity susceptibility, since Eq. (fJTJ is valid in 
the full range of coupling as long as T <§C Ar, we only need 
to verify that for N odd the singlet-triplet excitation energy 
SEst is, as expected from the same argument, a universal 
function of (Tk/Ar): <5T st /Ar = T(T k /Ar) with the lim- 
iting behaviors 



F(x) 



1. 



x > 1 



-l/ln(x), x < 1 . 



(24) 



The strong coupling behavior follows directly from the dis- 
cussion in Sec. IIII A I The scaling behavior of Fgx at weak- 
coupling can, on the other hand, be obtained from a perturba- 
tive renormalization group argument: in the perturbative ex- 
pression for SEst Eq. (T% , replace the coupling constant J 
by its renormalized value J e s at the scale Ar. Within the 
one-loop approximation, J c n = J/[l — J^\n(D/An)] yields 



SE; 



ST 



1 - Jln{D/A R ) 



(25) 



Substituting the one-loop expression for the Kondo tempera- 
ture, Tk — D exp(— 1 / J), now yields the second line of ( f24b . 
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FIG. 9: (Color online) Universal form of the singlet-triplet gap (N 
odd) in the absence of mesoscopic fluctuations ("Kondo in a box 
model"). Note the excellent data collapse for a wide range of bare 
parameters upon plotting the extracted gaps as a function of 7k/ Ar. 
Inset: Pqmc(P) for D/Ar — 41. Circles are data for fixed J = 
0.05, 0.10, 0.13, 0.15, 0.19, 0.23, 0.26, and 0.30 (right to left). 
Lines are one parameter fits to the two state singlet-triplet form used 
to extract the values of the gap. 



In the crossover regime, we find F(x) through continu- 
ous time quantum Monte Carlo (QMC) calculations using 
a modification of the algorithm presented in Ref.|26j with 
updates which maintain the number of particles (canonical 
ensemble). 63 To extract 6Est, we measure the fraction P of 
states with (S^ ot ) 2 — 1 visited in the Monte-Carlo sampling at 
temperature T. For a fixed J and large j3 — 1/T, P(j3) can 
be excellently fit to the form 2/(3 + e /3l5 - EsT ) valid for a two- 
level singlet-triplet system. Repeating this procedure yields 
5Est for a variety of J and Ar . 

Fig. [9] shows the results of our calculations, plotted as a 
function of Tk/Ar. We emphasize three features: (i) The 
inset shows that the fit of our QMC data to the simple two- 
level singlet-triplet form is indeed very good, (ii) The limiting 
behaviors of Eq. (l24i i are clearly seen, (iii) Data for a wide va- 
riety of bare parameters is shown; the excellent collapse onto 
a single curve in the main figure is a clear demonstration of 
the expected universality. 

We now turn to the local susceptibility and explore the ex- 
pected scaling ansatz 



f 



T T 

Ar"' Tk 



(26) 



for variety, we use the grand canonical ensemble. By fix- 
ing the chemical potential in the middle of the spectrum of 
the reservoir, particle-hole symmetry ensures that even in the 
presence of the Kondo coupling the mean number of particles 
(N) has a fixed parity: if p is aligned with a level, (N) is odd, 
while if /i falls exactly between two levels, (N) is even. 

Figs.[T0] and [TT] show our QMC results. In Fig.[l0l we 
demonstrate the expected scaling by showing the susceptibil- 



ity as a function of T/Tk for fixed TJ Ar [ie. "slices" of the 
function / in d26l ) are shown]. A wide variety of bare parame- 
ters are used in order to map out the full cross over, and a good 
data collapse is found. Note that in the low temperature limit, 
the universal function for (N) even is substantially different 
from that for (N) odd. 

One technical point concerning the determination of Ik: 
We expect that Tk is not affected by finite size effects (in the 
absence of mesoscopic fluctuations) - it is determined by the 
mean density p, the bandwidth D, and the coupling J in the 
same way as the "bulk problem". Indeed, this is explicitly ver- 
ified in Refs.l25landl26l However, for practical numerical pur- 
poses the 2-loop perturbative formula T^ e = D \f~Tp~ e~ x I J p 
works only approximately for J p > 0.20, and so we need to 
fit Tk numerically in this regime. Since the overall scale of 
the Kondo temperature is arbitrary, we set T5* = yP crt for 
Jp = 0.19. Having fixed this value, we can then obtain for 
all other couplings by choosing it to get the best data collapse 
in Fig. [lO] If we were truly in the perturbative regime, this 
extracted value would always agree well with the perturbative 
formula. The comparison between and J^ crt is made in 
the inset of Fig.[l0] As expected, for stronger couplings there 
are deviations from the perturbative value. After this fit, there 
are no further free parameters used to analyze the data. 

While the way the data is presented in Fig.[l0]is convenient 
for demonstrating the universal scaling Eq. ( f26b . it is more 
natural to present the data for fixed values of Tk/Ar since 
this corresponds to a fixed geometry or Hamiltonian. This is 
done in Fig.QT| The curves for odd and even (N) agree for 
T/Ar > 0.2. For lower temperature, the susceptibility satu- 
rates for (N) odd, while it shows the expected Curie law for 
(N) even. When scaled by Tk [panel (b)], the curves for both 
parities follow the bulk universal curve as temperature de- 
creases until they separate from each other when T j Ar « 0.2. 



D. Limiting regimes and interpolations for T <C Ar 

To close this section, we come back to the limiting behav- 
iors of the susceptibility discussed earlier. We show in detail 
how the interpolation between them takes place, as well as 
the relation between the results for the canonical and grand- 
canonical ensembles. Fig. [T2l shows our QMC results yet an- 
other way: for fixed T/Ar = 0.1 but over the whole transi- 
tion regime of Tk/Ar and for both ensembles. As expected, 
the results from the canonical ensemble also satisfy the scal- 
ing ansatz. We see moreover that the weak and strong cou- 
pling limits of the canonical QMC data are reasonable: it is 
in good agreement with the expressions Eq. d22"i i for the weak 
coupling regime, and tends to the bulk zero temperature limit 
/ x (0) for large Tk/Ar. (Note, however, that the numerics 
were not pushed as deep into the limiting regimes as for <5Tst 
in Fig.[9]). 

The strong dependence of the cross-over behavior on the 
parity is unexpected. The odd case has a featureless transi- 
tion and reaches the large Tk/Ar limit from below. In con- 
trast, the susceptibility in the even case overshoots the large 
Tk/Ar value and then approaches / x (0) from above. As the 
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FIG. 10: (Color online) Rescaled susceptibility Xuniv = ^loJxioc [in units of (q/j-b) 2 /4] in the grand canonical ensemble as a function 
of T/Tk while keeping the ratio T/Ar fixed. The four panels correspond to T/Ar = 0.5, 0.2, 0.1, and 0.05. In each panel, the light 
(yellow online) and dark symbols correspond respectively to odd and even mean number of particles in the bath. Note the collapse of the 
data on universal curves, and the substantial difference between the even odd cases at low temperature. Inset: Comparison between the Kondo 
temperatures extracted from our fits and the two-loop perturbative formula. As expected, the agreement is excellent for J <C 1. At larger values 
of J, the two-loop formula tends to underestimate T K . In each panel there is data for £>/A R = 20(21), 40(41), 80(81), 160(161), 320(321) 
for the odd (even) case, and for each case we have used 7 = 0.19, 0.23, 0.26, and 0.30. 




FIG. 11: (Color online) Rescaled susceptibility Xuniv in the grand canonical ensemble. Same data as in Fig.[l0] but now plotted for fixed 
values of Tk/Ar, i.e. a fixed Hamiltonian. Points connected by solid (dotted) lines have 161 (160) sites in the bath, resulting in an odd (even) 
value for (N). Results are plotted in units of (a) Ar and (b) Tk. Magnitude of coupling: from top to bottom in (a) and outside to inside in (b), 
7 = 0.19, 0.23, 0.26, and 0.30; these correspond to Tk/Ar = 0.18, 0.54, 1.03 and 2.02 (Tic used here determined as in Fig.fTot. 



even and odd canonical results are so different, this provides 
an excellent test of the connection between the canonical and 
grand-canonical results discussed in Section [VBl the grand- 
canonical result in the even case agrees with the average of 
the two canonical results (green line). As a consequence, the 
grand-canonical susceptibility for even (N) differs consider- 
ably from the even iV canonical result even for fairly large 
Tk/Ar. 

For weak coupling, although the two expressions in 
Eq. d22b appear quite different, they lead to similar numerical 



values for a very large range of T/Ar, and in particular for 
the value 0.1 used in Fig. [12] As a consequence, the canonical 
and grand-canonical data in the odd case do not differ very 
much in the weak coupling regime. Note however that this is 
somewhat emphasized here by the fact that Tkxioc is plotted 
rather than simply the susceptibility itself. 
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FIG. 12: (Color online) Comparison of canonical and grand-canonical ensemble results for the rescaled susceptibility Xuniv = ZiocXioc [in 
units of (ij/is) 2 /4] at the fixed temperature T/Ar = 0.1 for an (a) even and (b) odd number of particles. Symbols: Quantum Monte Carlo 
results for the canonical (circles) and grand-canonical (triangles) ensembles. Solid line: Bulk zero temperature limit. Dashed lines: Weak 
coupling expression Eq. d22b : for odd case, uses QMC values of SEst- Doted line: Weak coupling expression Eq. d22b using the asymptotic 
expression J24b for SEst- Dash-dotted lines: Average of the even and odd canonical results to compare with the grand-canonical result [using 
the numerical quantum Monte Carlo data in (a) but the asymptotic expressions Eq. d22t in (b)]. 



VI. CONCLUSIONS 

Our focus in this paper has been on the many-body spec- 
trum of a finite size Kondo system. We have in mind a "mag- 
netic impurity" - either a real one or an effective one formed 
by a small quantum dot - coupled to a finite fermionic reser- 
voir (Fig.Q]). Such a system can certainly be made with current 
technology; indeed, using quantum dots, a tunable connection 
between the small dot (S) and reservoir (R) could be made, al- 
lowing a direct investigation of the parametric evolution of 
properties as function of the coupling. The emphasis through- 
out the paper is on experimentally observable consequences. 

We start with a theorem for the ground state spin of the 
combined R-S system. Because a crossing of the ground state 
is forbidden, this can be obtained from simple perturbation 
theory - the result for different cases is given in Table I. The 
theorem is a straightforward extension of the classic theorems 
of Mattis 43 and Marshall. 44 

A schematic picture of the spectrum of low-lying states as a 
function of the coupling J can be constructed using perturba- 
tion theory and plausibility arguments. Fig s.[5] and [6] show re- 
sults for the screened and under-screened cases, respectively. 
In this respect we are greatly aided by having a perturbation 
theory available not only at weak coupling but also at strong 
coupling (Nozieres' Fermi liquid theory). 

The first observable property that we focus on is the nonlin- 
ear I-V curve of such an R-S system. Using a rate equation 
approach, we find the differential conductance as a function of 
the bias voltage for a number of cases (Figs.|7]and[8]l, with cer- 
tain simplifying assumptions so that the identification of the 
different transitions is clear. The key result is that the splitting 



with magnetic field combined with the magnitude of dl/dV 
can be used to deduce the spin of the low-lying excited states. 
Thus an experiment could obtain the information needed to 
compare with the theoretical schematic picture. 

The second observable that we treat is the magnetic sus- 
ceptibility; in the low temperature limit, this is simply related 
to the low-lying states of the system. We study the general 
behavior of both the impurity and local susceptibility, finding 
that they are markedly different. An extensive example illus- 
trates the general features: A quantum Monte Carlo calcula- 
tion yields results for the singlet-triplet energy gap (impurity 
susceptibility) and the local susceptibility in the "clean Kondo 
box" model in which the levels are equally spaced and all lev- 
els couple to the impurity with the same amplitude. Results 
in Figs.l9lfl~2lshow the expected universality in this model, the 
strong even-odd effects, and the difference between using the 
canonical and grand-canonical ensemble. 

Clearly the finite size spectrum of engineered many-body 
systems is a rich area for future experiments. We hope that 
our schematic arguments plus results for two observable prop- 
erties will persuade researchers to undertake them. 
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APPENDIX A: MARSHALL'S THEOREM 

Theorem: If a Hermitian matrix H and a basis set \<f> a ) sat- 
isfy two conditions, 

(i) {4> a \H\<j)p) < for a ^ 13, and 

(ii) for every pair a and (3, (4> a \H n \(t)p) ^ for at least one 
integer n, 

then the ground state of H is unique (i.e. non-degenerate). 

We prove the theorem in a few steps: 

Step I: First, we prove that if the conditions above are satis- 
fied, then the following simple result is true: if J2 a /a|0a) is 
a ground state, then so is ^ Q \f a \ \(f> a ). To show this, consider 
the following "Perron-Frobenius" inequality: 

\M(^\H\h) < E fMtMfy) ■ (Al) 

a,0 a, (3 

Note that the two sides of the inequality differ only for terms 
with a 7^ (3. All such terms on the left-hand side are negative 
because of condition (i) above, while on the right-hand side 
the sign may be positive, depending on the sign of f a ; hence, 
the inequality holds. Physically, the right and left side of the 
inequality are, in fact, the expectation values of the energy 
in the ground state ^ a fa\4>a) and in J2 a l/a I 10a), respec- 
tively. Since the expectation value of the Hamiltonian in the 
ground state is always the lowest, D Q \f a \ \<f> a ) must also be 
a ground state. 

Step II: We now show that any ground state described by f a 
must have f a > 0. For convenience we separate H = Htj + 
Hqd, where Hu and Hqd are the diagonal and off-diagonal 



parts of H with respect to the basis \(f> a ). In this language the 
Schrodinger equation expressed for the states ^ Q f a \4> a ) and 
La l/a I l^a) becomes 

(0a|ff D |0a)/a - ^ I (4><* l#OD 1 |//3 = E G f a (A2) 


(<t> a \ H n\<l>c*)\f*\ - K<M#od|0/3>| \fp\ = E G \f a \ (A3) 



where Eq is the ground state energy. These two equations 
combine to give 

|«0a|tfD|0a)-£ G )/a| = | (<£a |flb 1 4> a ) ~ E G | | f a \ . 

(A4) 

Note, however, that (4> a \H-£>\(f>a) ~ Eq > 0, because \<p a ) 
cannot have a lower energy expectation value than the ground 
state. Thus, f a >0 follows. 

Step III: The next step is to use condition (ii) from the the- 
orem to show that the stricter condition f a >0 holds. This is 
most easily seen in Eq. (IA2b : if f ai — for one a±, then, using 
condition (ii), all the f a must be zero. Since this cannot be the 
case, all f a ^ 0. Thus, combining with the result of Step II, 
we conclude that f a > 0. 

Step IV: We have thus shown that every ground state of 
H has a positive definite expansion in the basis set \4> a ). 
Since there cannot be two simultaneously orthogonal and pos- 
itive definite vectors, the ground state of H must be non- 
degenerate. This proves the theorem. 
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